Effect of Dynamical Coulomb Correlations on the Fermi Surface of Nao.3Co02 
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The t2g quasi-particle spectra of Nao.3Co02 are calculated within the dynamical mean field 
theory. It is shown that as a result of dynamical Coulomb correlations charge is transfered from 
the nearly filled Cg' subbands to the aig band, thereby reducing orbital polarization among Co t2g 
states. Dynamical correlations therefore stabilize the small Cgi Fermi surface pockets, in contrast 
to angle-resolved photoemission data, which do not reveal these pockets. 
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The intercalated layer compound Na2;Co02 has been 
intensely studied during the recent years since as a func- 
tion of Na doping it exhibits a variety of fascinating 
properties. The relevant valence bands near the Fermi 
level consist of Co t2g states with occupancy 3d^+^ . For 
X « 0.50 — 0.75 an unusually large thermopower is ob- 
served 0. For X « 0.3 hydration gives rise to a supercon- 
ducting transition at 4.5 K Q. In a narrow region near 
a; = 0.5 the material undergoes a metal insulator transi- 
tion 01 . The end-member at a; = with a single hole 
per Co atom is believed to be a Mott insulator. At x = 1, 
since the filled t2g bands are separated by about 1.5 eV 
from the empty Cg bands, one finds a band insulator. 

Despite considerable experimental and theoretical ef- 
fort, fundamental electronic properties of Naa:Co02 such 
as the qualitative topology of the Fermi surface remain 
controversial and are not well understood. In the metallic 
phase the hexagonal Fermi surface predicted by band the- 
ory within the local density approximation (LDA) con- 
sists of a large hole pocket centered at T and six small 
hole pockets near K [5j (see Fig. 1). Because of the lay- 
ered structure the t2g levels split into a aig level and dou- 
bly degenerate eg' levels. The large Fermi surface stems 
mainly from the aig bands while the small hole pock- 
ets have predominantly eg' character. The existence of 
these hole pockets is thought to be crucial for the under- 
standing of the superconducting phase |^ . On the other 
hand, angle resolved photoemission spectra (ARPES) for 
X = 0.7 and x = 0.6 Q provide evidence only for aig 
bands crossing the Fermi level. These data suggest that 
the egi bands are filled due to inter-orbital charge trans- 
fer not described within the LDA. New ARPES data for 
X — 0.3 • • • 0.7 and x — 0.3 (1^] also do not show the 
small eg' hole pockets. 

To address the question of possible modifications of 
the Fermi surface via Coulomb correlations not included 
in the LDA, several theoretical studies were carried out 
for the half- metallic, magnetic region near x — 0.7 [l]| . 
using the LDA-I-U approach For a relatively large 
local Coulomb energy ([/ > 3 eV), the eg' subbands are 
indeed filled. In the metallic phase, however, in par- 
ticular, for Coulomb energies comparable to the single- 



electron band width, it is well known that correlations 
have an important dynamical component not captured 
within the LDA-I-U. 

The aim of this work is to elucidate the possibility of 
modifying the Fermi surface of Nao.3Co02 via dynamical 
Coulomb correlations. More specifically we focus on the 
charge transfer between the paramagnetic t2g subbands. 
Since only the total electron number is conserved when 
correlations are taken into account, the occupations of in- 
dividual subbands may vary with the strength of the local 
intra- and inter-orbital Coulomb energies. In a single- 
band picture, the key effect of dynamical fluctuations is 
the spectral weight transfer from the quasi-particle peak 
near Ep to the incoherent satellites associated with the 
lower and upper Hubbard bands. In a multi-band ma- 




FIG. 1: Upper panel: Fermi surface of Nao.3Co02. Lower 
panel: tight-binding fit to LDA bands; Ef ~ 0. Solid curves: 
even bands; dashed curve: odd band. The Cg' states above 
Ef give rise to the small hole pockets of the Fermi surface. 
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terial, this spectral weight transfer is orbital dependent, 
opening the possibility of redistributing electronic charge 
among the valence orbitals and modifying the shape of 
the Fermi surface. 

To investigate these multi-band correlation effects we 
use the dynamical mean field theory (DMFT) combined 
with the Quantum Monte Carlo (QMC) method 
The remarkable result of this work is that in the metallic 
domain of Na2,Co02 near x = 0.3 dynamical correlations 
shift charge from the Cgi states to the aig bands, thereby 
stabilizing the eg> hole pockets and slightly reducing the 
aig Fermi surface. The overall topology of the Fermi 
surface remains the same as in the LDA. 

Fig. 1 shows a tight-binding fit to the Nao.3Co02 bands 
calculated within the LDA and using the linearized aug- 
mented plane wave (LAPW) method. We consider the 
paramagnetic phase observed experimentally at a; = 0.3. 
Since we are interested in the qualitative issue of charge 
transfer between t2g subbands the weak dispersion along 
the c-axis is neglected. The full band structure involving 
Co 3d and O 2p states is down-folded to a 3 x 3 Co t2g 
tight-binding Hamiltonian in which on-site energies and 
hopping integrals represent effective energies accounting 
for direct Co-Co and indirect Co-O-Co interactions. In- 
cluding three neighbor shells, with dda, ddir and ddS ma- 
trix elements, an excellent fit to the LAPW band struc- 
ture is achieved. The details of the tight binding model 
and down- folding procedure will be given elsewhere |l4l |. 

Because of the planar structure of the system, it is 



convenient to transform the dx 



orbitals into ai„ and 



eg' states, where aig = {dxy+dxz+dyz)/V^, egi = {d^z- 
dyz)/V^, eg2 = i2dxy - dxz ~ dyz)/V^- Quantities such 
as the local density of states and local quasi-particle self- 
energy are diagonal in this representation. For instance, 
the aig and eg' density of states are pa = pa + "^Pij, 
Pe = Pii-Pij, where p^i (pn = p22 = P33) and pi^j {pi2 = 
pi3 = P23) are the diagonal and off-diagonal elements of 
the t2g density of states matrix. 

Fig. 2 shows the aig and eg' quasi-particle spectra 
for Nao.3Co02 as calculated within the DMFT. The lo- 
cal Coulomb interaction defining the quantum impurity 
problem is characterized by intra- and inter-orbital ma- 
trix elements U ~ Uu and U' = Ui^j = U — 2 J, where 
J is the Hund's rule exchange integral. The value of U 
for the entire ^23, eg manifold (total width about 4 eV) 
was estimated at 3.7 eV [31. For the narrower t2g bands 
a smaller value should be more appropriate to account 
for screening involving empty Cg states. Since accurate 
values of U and J are not available, DMFT calculations 
for several values were carried out to study their effect 
on the inter-orbital charge transfer. The temperatures 
were T « 385 and 770 K corresponding to about 30 and 
60 meV thermal broadening. Up to 10^ sweeps were done 
in the QMC calculations. The quasi-particle spectra are 
obtained via maximum entropy reconstruction . 

The quasi-particle spectra show the characteristic band 
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FIG. 2: Quasi-particle spectra for metallic Nao.3Co02 cal- 
culated within the DMFT (solid curves) for (7 = 3 eV, 
J = 0.8 eV; T = 770 K, Ef = 0. Dashed curves: LDA 
local density of states, (a) a\g states; (b) Cgi states. 



narrowing near Ep caused by dynamical correlations and 
the transfer of weight from the coherent to the incoherent 
spectral region. In the slightly narrower eg' band corre- 
lations are strong enough to give rise to a lower Hubbard 
band. Also noticeable is the substantial lifetime broaden- 
ing of valence states due to creation of electron hole pairs. 
The occupations of these distributions are: Ua — 0.853, 
Uf, — 0.899, which should be compared to the LDA values 
Ua = 0.797, Ue = 0.927. 

To illustrate the variation of the t2g subband occu- 
pations with local Coulomb and exchange energies, we 
show in Fig. 3(a) the trend obtained within the DMFT 
for U' = C//2, J = C//4 and [/' = [7, J = 0. In both 
cases charge transfer proceeds from eg' to aig, i.e., orbital 
polarization is reduced. We have evaluated the DMFT 
quasi-particle spectra both in the non-diagonal t2g and 
diagonal aig, eg' representations. Both versions are in 
excellent agreement, indicating that in the present sys- 
tem non-diagonal coupling among t2g states is fully taken 
into account within the aig, eg' representation. 

The variation of the oig, eg' subband occupations with 
U derived within the DMFT differs from the one found 
in the LDA+U, as illustrated in Fig. 3(b). While dy- 
namical correlations lead to reduced orbital polariza- 
tion for J = U/4 and J = 0, the LDA-I-U treatment 
gives this trend only for J = ?7/4; J = yields the 
opposite effect. This dependence of the orbital polar- 
ization on the ratio J/U within the LDA-I-U follows 
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FIG. 3: Occupations of aig and Cgi subbands of Nao.3Co02 
as a function of U. (a) DMFT; (b) LDA+U. Solid dots: J = 
17/4, empty dots: J = Q. The lines are guides for the eye. 



from the Hartree Fock self-energy. For a paramagnetic 
t2g complex with one-fold aig and two-fold eg/ sub- 
bands, the orbital dependent potential is given by [l^ : 



LDA-hU 



= Hff = S{U - 5 J), with S = {ue- na)/3. 
The diagonal term = 5n{U — 2J) gives an over- 

all energy shift, where n = {ua + 2ne)/3. Within the 
aig, Cgi basis, the self-energies are Sjj'^ = -I- 2S^-^, 
E^^ = S^^ - S^F. Subtracting the diagonal term the 
shifted band energies are: f-'a{k) = ea{k) + 26{U — 5 J) 
and ee(fc) — (-e{k) — 5{U — 5J). Thus, for rie > Ua 
the aig (cg') energies are shifted up (down) as long as 
J < U/5. For J > U/5 this trend is reversed. For real- 
istic U = 3.0 • • • 3.5 eV, J « 0.8 eV, the latter condition 
is satisfied, implying reduced orbital polarization also in 
the static limit. 

Comparing the static self-energy with the one obtained 
in second-order perturbation theory we find that, be- 
cause of the small difference between and rig, for in- 
creasing U static correlations are rapidly dominated by 
dynamical correlations. Also, the more compact eg' den- 

(2) (2) 

sity of states ensures that Re Se (w = 0) > Re Tia {lj — 
0) for J = and J — U/A, implying diminishing orbital 
polarization [l^ . Thus, static and dynamical correla- 
tions exhibit qualitatively different dependencies on the 
ratio J/U. In addition, the reduced orbital polarizations 
found for J — L//4 in the static and dynamical cases arise 
for different physical reasons. 




FIG. 4: Occupations of 3-band model: elliptical density of 
states; band fillings as in Nao.3Co02; DMFT. Solid dots: J = 
?7/4, empty dots: J — 0. The lines are guides for the eye. 



The results in Fig. 3 demonstrate that in the metallic 
phase it is crucial to include dynamical correlations. Ev- 
idently, the new degrees of freedom generated by quasi- 
particle interactions, such as spectral weight transfer be- 
tween low and high frequencies, relaxation shifts (band 
narrowing) and decay processes, which are beyond the 
LDA and LDA+U treatments, contribute to the charge 
balance between non-equivalent orbitals. It is to be ex- 
pected that these dynamical correlations also depend on 
the single-particle bands, i.e., on the available density 
of occupied and unoccupied states involved in excitation 
processes. To illustrate this point we show in Fig. 4 the 
subband occupations for an analogous 3-band model with 
elliptical density of states and total occupancy 5.3. Ad- 
justing the chemical shift to simulate the aig, Cg' occu- 
pations of Nao.3Co02 we find increasing orbital polariza- 
tion for J = but weakly varying rii for J = U/A [l^. 
Clearly, the inter-orbital charge transfer depends on the 
combination of single- and many-particle interactions. 

The correlation induced band narrowing seen in Fig. 2 
is consistent with infrared optical data jj| which show 
transitions within the t^n manifold at lower energies than 
predicted by the LDA Substantial band narrowing 
is also observed in all ARPES measurements 0jS0E|- 
On the other hand, the present results are in conflict 
with the ARPES data as far as the filling of the Cg' 
bands is concerned. So far, we have no explanantion for 
this qualitative discrepancy between our state-of-the-art 
DMFT results and the ARPES data. The interpretation 
of photoemission spectra can be complicated due to sur- 
face induced single-particle and correlation features [l^l ■ 
Moreover, matrix element effects could play a role. For 
instance, along FK (chosen as x axis) the aig and eg2 
bands are even with respect to zL{y,z) where z defines 
the surface normal (see Fig. 1). In contrast, Cgi is odd. 
These states therefore couple differently to the incident 
photon field. To identify the orbital character of the t2g 
bands near Ep the use of various polarizations is recom- 
mended. 
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It is interesting to compare the influence of dynamical 
correlations on the subband occupations of Nao.3Co02 
to those found in other multi-band materials. The 
valence bands of the layer perovskite Sr2Ru04 also have 
aig, eg' symmetry, but are 2/3 filled. Because of the 
smaller width of the eg' bands DMFT predicts a charge 
transfer to the aig band 0|, shifting the ojo van Hove 
singularity closer to Ep than in the LDA |21|. The ba- 
sic shape of the Fermi surface is preserved. The trans- 
fer from Cg' to aig agrees with the one found here for 
Nao.3Co02. A similar inter-orbital charge transfer was 
recently obtained in DMFT calculations for BaVSa [2^ . 
In the orthorhombic phase this material exhibits a 
wide aig band and weakly occupied, narrow eg' bands. 
Dynamical correlations were found to cause a transfer of 
electrons from aig to eg' states, thereby reducing orbital 
polarization. As in the present system, the same trend 
was obtained for small and large Hund's rule exchange 
terms (J = U/7, U/A). On the other hand, as recently 
shown by Pavarini et al. [l^ for several 3d^ perovskites, 
non-diagonal t2g coupling caused by octahedral distor- 
tions decreases orbital fluctuations and enhances insulat- 
ing behavior. Also, Manini et al. studied the Mott 
transition in a model consisting of two equal subbands 
with unit total occupancy. With increasing chemical 
shift, causing more diverse band fillings, the tendency to- 
wards a Mott insulator was found to be enhanced. Since 
in their case J = 0, these results are consistent with the 
enhanced orbital polarization shown in Fig. 4 for J = 0. 
These various trends underscore the subtle nature of cor- 
relations in multi-band systems, and their remarkable 
sensitivity to various system parameters, in particular, 
the shape of the density of states and the size of J/U. 

In summary, we have explored the modification of the 
Fermi surface of Nao.3Co02 as a result of dynamical 
correlations. The main effect of these fluctuations is 
the spectral weight transfer from the quasi-particle peak 
near Ep to the incoherent part of the spectrum. In a 
multi-band material, these dynamical processes are or- 
bital dependent. In particular, since they depend on or- 
bital occupations, they can give rise to charge transfer 
between subbands induced via hybridization and inter- 
band Coulomb interactions. The highlight of the LDA 
band structure of Nao.3Co02 is the distinctly different 
filling of the t2g bands, with the nearly full eg' subbands 
yielding characteristic hole pockets of the Fermi surface. 
Accounting for correlations within the DMFT we have 
shown that electronic charge is shifted from the eg' sub- 
bands to the aig band, thus slightly enlarging the small 
eg' hole pockets and reducing the main aig pocket cen- 
tered at r. Further studies are needed to reconcile these 
theoretical findings with the photoemission data. 
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